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Numerical analysis of the spectrum of large finite size Luttinger liquids (g < 1) in the presence 
of a single source of 2k f scattering has been made possible thanks to an effective integration of high 
degrees of freedom. Presence of irrelevant operators and their manifestation in transport are issues 
treated independently. We confirm the existence of two irrelevant operators: particle hopping and 
charge oscillations, with regions of dominance separated by g = 1/2. Temperature dependence of 
conductance is shown to be controlled by hopping alone, giving G ~ j^/s- 2 Frecuency dependence 
is affected by both operators, giving G ~ lo 2 ^ 9 ^ 2 for g > 1/2, and G ~ uj 2 for g < 1/2. 

I/-) ' PACS numbers: 72.10.Bg, 73.20.Dx, 73.40. Gk 

0^ ' The search for metallic non- Fermi liquid behavior, motivated by normal state properties of high T c superconductors 
jlj, has renewed the interest in one-dimensional (Id) interacting system. Interactions endow particles in Id with 
OJQ' properties that are qualitatively different from the quasi-particle Fermi liquid picture, requiring a new description 
under the name of Luttinger liquids (LL) g|. Transport is severely affected and, for instance, the conductance of a 
' perfect channel passes from being universal to depend on the interaction strength ||. More striking is the effect of 
, even a single impurity: depending on the host conductance, the imperfection can either transmit or reflect perfectly 
Although there have been great advances in the nanofabrication of Id structures Q, the realization that edge 
states of quantum Hall effect liquids are chiral LL |jj has opened new perspectives to this field. Now it seems possible 
to give experimental reality [|| to idealizations such as a single source of 2k f scaterring in an otherwise perfect LL. 
A LL in the presence of a localized source of 2k f scattering is described by 
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j where 9 and II are canonical conjugate fields related to charge and current respectively, g is the conductance in the 
absence of perturbation, and the impurity scattering is given by the last cosine term. For g < 1 (repulsive case) the 
perturbation is relevant MM, effectively splitting the system into two parts, with the phase at the impurity location, 
9 = 9{x = 0), pinned in any of the equivalent minima of the cosine term. The pinning of 9 blocks transport 
driving the conductance to zero B . At low energy, 9 owes its residual dynamics to hopping processes between cosine 
minima (particle transport across the junction). This hopping is an irrelevant perturbation with scaling dimension 
l/<7, giving the current-current correlator at the junction the imaginary time (r) dependence t~ 2 ^ 9 . Standard scaling 
for Kubo formula would then provide the following temperature (T) or frecuency (u>) dependence to the conductance: 
G ~ r/ 2 /9- 2 , with r) = uj or T, as Kane and Fisher (KF) have shown H, 

Based on the Self- Consistent Harmonic Approximation (SCHA) f8|j, which substitutes the effect of the cosine in 
Eq. (Q) by a quadratic term a 9 2 , Guinea et al. fTo)] have qualified this picture. They claim that, in addition to 
the hopping between minima, the SCHA suggests the existence of another irrelevant perturbation of dimension 2, 
describing oscillations of o around its pinned value and providing the junction with a capacitance. (A density-density 
coupling between both sides of the already split system would be an equivalent description.) This new irrelevant 
operator (dominant for g < 1/2) would give rise to a ^-independent r -4 contribution to the junction current-current 
correlation. Scaling then dictates a io 2 dependence for the conductance, replacing the KF result uj 2 / g ~ 2 when g < 1/2. 

The case for this new operator seems strong. Based on a perturbative analysis, Chamon et al. [ pT[ have found 
evidence for the t~ 4 dependence of the current-current correlator for g < 1/2. The same authors |12| have been 
forced to invoke the presence of a density-density coupling to understand features of the noise spectrum, and Tsvelik 
fyl has deduced the existence of a capacitance from Bethe ansatz techniques. On the other hand, the same scaling that 
gives the uj 2 law for the dynamic conductance would predict a T 2 temperature dependence for the static conductance 
JT(| replacing the KF prediction for g < 1/2. Unfortunately this last result contradicts the exact solution found for 
g = 1/3 from Bethe ansatz techniques fL4| , which is in agreement with the KF analysis G ~ T 2 I 9 ~ 2 . This confronts 
us with a puzzling situation. 

Naive scaling is not guaranteed to work, and to settle the problem described above one would like to separate the 
following two issues: (i) presence of irrelevant operators (nature, dimension, and region of dominance), and (ii) their 
manifestation in physical properties. With this in mind, we have adopted a finite size approach in this paper. The 
idea is well known: solving H in a finite size L amounts to probing the system with a energy scale set by 1/L. The 
low energy spectrum of the (scaled) Hamiltonian LH becomes scale invariant in the limit L = oo (fixed point). The 
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presence of irrelevant operators and their scaling dimensions can be read off from the L depedence of the approach 
to the fixed point. 

Unfortunately, this simple idea faces severe implementation problems. The number of degrees of freedom grows 
gigantically with size, and there is no real chance of solving exactly a finite system of even modest size. One should 
resort to approximations that keep a manageable number of states while ensuring that they truly describe the low 
energy sector of interest. In our case, the use of the SCHA as an intermediate step will help us handle this problem 
in a transparent and reliable way, as we will see. 

It is well known that, under the guise of a variational approach, the SCHA performs a perturbative renormalization 
p[-|lO[: replacement of the cosine term in Eq. (|l|) by a self-consistent a9\ provides the following dependence: a ~ 
vr~~ 9 \ This behavior is precisely what would have been obtained for the running coupling u upon integration 
of its renormalization group (RG) flow up to the point where u equals the cut-off energy, signalling the end of the 
perturbative regime. Thus, the low energy states of the SCHA have a built-in effect of the high energy modes. 
Although the SCHA flows to the same fixed point expected for the original problem ( 9 a is also pinned), it loses 
memory of the discrete (global) translational invariance (9^9 + mr) of H . It is clear that the SCHA can pin 9 Q 
at any of the translation-equivalent minima, suggesting the following procedure: replicate the SCHA states for each 
well, and use these states as a basis for the diagonalization of the original Hamiltonian. 

The steps in the precise implementation of this idea are the following, (i) The Hamiltonian H is regularized in 
a linear chain of L sites with periodic boundary conditions, the cosine perturbation affecting only one lattice site, 
(ii) The SCHA is performed for every L, providing a basis which is replicated at every equivalent well, (iii) The 
regularized hamiltonian is scaled, h(L) = HL/n, and diagonalized in the replicated SCHA basis, keeping a fixed 
number of states (~ 100) per well and taking advantage of the 9 translational symmetry. 

It is important to stress that the true Hamiltonian is diagonalized in a restricted basis provided by he SCHA as an 
intermediate step. The effective integration of high modes implied by the SCHA will prove crucial for the success of 
this approach. Matrix elements and overlaps pose no severe computational problem owing to the harmonic nature of 
the basis. The calculation has all the ingredients of tight-binding band structure problem, and the spectrum adopts 
the form of enery (e) versus crystalline momentum (q), associated with the (discrete) ^-periodicity. 9 being charge, q 
can be identified with flux piercing the ring, enabling transport properties to be read off from e(q) p5[ . 

We have performed calculations following the above described procedure for values of g in the interval 1/3 < g < 1, 
with lattice sizes in the range 100 < L < 10 7 , and values of u bounded by u < 10 _1 (much less than the energy 
cut-off). The number of states per 9 well is 100, having checked that doubling this number produces no significant 
changes. The first important point to remark is the (it, L) dependence of the spectrum of the scaled h(L) through the 
combination uVC 1- ff)i s allowing us to measure the effect of the perturbation by the following scaling variable: 

u(u,L) = u{2/L) g - 1 (2) 

Trading 1/L for temperature, this is the same scaling variable found by KF || to describe their conductance results. 
This scaling is a consecuence of the RG flow for the relevant coupling u. In our scheme, this property comes through 
the mentioned dependence a ~ ^/(i-a) Q f the SCHA, emphasizing again its implicit integration of high degrees of 
freedom. 

Representative results for the spectrum are shown in Fig. 1, where the two lower bands £1,2(9) of h(L) for g = 0.4 
are plotted for increasing values of u. (Similar results are obtained for all values of g.) The spectrum evolves smoothly 
from the fixed point u = + to that of u — > 00, where (flat band) energies and degeneracies are, of course, those 
of a broken chain. Although our interest is the approach to this last fixed point, the analysis of the departure point 
u ~ provides a critical test of our treatment. In the left panel of Fig. 1 we have superimposed the exact result for the 
(folded) parabola of current-carrying states of a perfect LL with our results. The agreement is perfect, and one should 
emphasize that the scheme of our calculation is the same for all values of u. That means, in the solid-state language, 
that we are reproducing the free particle dispersion starting out of a tight- binding basis. Following this parabola 
upwards in energy, discrepancies between exact and calculated results appear, due to the limited basis employed (this 
happens around e ~ 5 in the present case). This provides us with a quantitative measure of the validity of our results. 
All calculations presented here are for energies within this (u — 0) confidence limit, even though it is clear that this 
energy window will widens with increasing u. A further point to notice in the small u limit is that the Bragg gap, 
Seb = £2(71") — £1,(71"), opens with a dependence Sen ~ u ~ L 1 ^ 9 . This is the expected perturbative behavior for the 
2kf perturbation of a LL, which grows with scaling dimension g. This agreement in the nominally worse situation for 
a scheme based on a tight-binding approach, gives us confidence in the quality of our approach. 

Now we study irrelevant operators. The most conspicuous effect in the approach to the u = 00 fixed point is the 
lack of band dispersion. This reflects the vanishing particle hopping across the junction (tunnel between 9 Q minima). 
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We characterize this residual hopping by the dispersion of the lowest band, and present its L dependence in Fig. 2 
(left panel). For all values of g it shows the dependence 

J Cl = ( Cl ( 7r )- ei (0))~ L 1 - 1 ^, (3) 

corresponding to a scaling dimension 1/g, in agreement with the KF analysis for the hopping across the junction. 

Let us now consider other clear feature of the spectrum in Fig. 1: the approach of the second band to its asymptotic 
value at the Brillouin zone center, 5e2 = (1 — 62(0)). Its L dependence is plotted in Fig. 2 (right panel) and summarized 
as follows: 

^ = 11-^(0)1-1^ W<i/S) ( 4 ) 

If tunnel between wells were the only irrelevant operator, then we would always have the first result of Eq. (|4|) This 
is certainly what happens for g > 1/2, while for g < 1/2 the behavior obtained is precisely that expected from the 
presence of an operator with scaling dimension 2 describing fluctuations of the pinned phase 9 , as explained before. 
These results confirm the picture for two irrelevant operators described at the beginning, with the value g = 1/2 
marking the boundary between regions of dominance. 

Now we study the consequences for transport properties, beginning with the temperature dependence of the static 
conductance. Although we cannot calculate the conductance of the infinite system as a function of T, we can 
circumvent this problem by defining a running temperature for each size, T — 2/L, so that the energy scale for a 
given size is also the temperature. We define the conductance for each size-temperature according to 

G(u,L 1 T = 2/L) = ir 2 <d 2 q t> T , (5) 

where its well known expression as a thermal average (<>t) of the second derivative of energy versus flux has been 
used ||). 

Although the G so calculated is not the true (infinite size) conductance as a function of T, it is clear that it will 
be a universal function of the scaling variable u [Eq. (0)], with the same scaling dependences. Results for G(u) are 
plotted in Fig. 3 (left panel) for g = 0.4. Notice the correct limit G(u — > 0) — > g, and a smooth crossover to the 
asymptotic behavior G{u — > oo ) -> u- 2 /9. This limit implies (remenber T = 2/L) that G ~ T 2 /^ 2 , in agreement 
with the KF result. Similar results have been obtained for all values of g, at both sides of g = 1/2. This should not 
surprise in view of Eq. (j^), where only the operator responsible for band dispersion can appear. 

Now we study the dynamic conductance, that is, the uj dependence of the (real) part of G (zero temperature). 
Placing a time dependent voltage (represented by a vector potential) right at the junction adds a term to the H 
proportional to the current density n o . Calculating matrix clement of n o , we can extract the power loss associated 
with a transition frecuency, whereupon G(w) can be obtained. We define a (scaled) oscillator strength j 12 (<z) for the 
vertical transition between states of the two lowest bands as 

/12(g)- ( 5 /2)| < g,2|n |<z,l >\ 2 L 2 , (6) 

and analyze its L dependence for very large sizes. The results are presented in Fig. 3 (right panel), where the 
(Brillouin-zone averaged) oscillator strength, f 12 =< /12(g) >q, is plotted for several values of g. The following 
asymptotic behaviors are obtained: 

f < f M > ~ / L ^ 2 ' 9 (1/2 <g<1) (7) 
fu =< /12(g) > g ~ j L _ 2 {g < 1/2) (7) 

It is a simple exercise to show that /12 ~ L~ v implies G(ui) — uf 1 . Thus, our results show that the dynamic 
conductance changes from G{uj) — lj 2 / 9 ~ 2 for 1/2 < g < 1, to G(u>) — lo 2 for g < 1/2. 

This is in agreement with the picture presented in the introduction, meaning that both irrelevant operators con- 
tribute to G(ui). For g > 1/2 hopping dominates and one obtaines the KF result. For g < 1/2 oscillations of the 
pinned O dominate, and G(uS) reflects the capacitance interpretation of this irrelevant peturbation. It is interesting to 
notice that a simplified analytical treatment that keeps only the two lowest states per well, shows that the oscillator 
strength comes from two contributions. One performs the transition for the states of the same well (no particle 
tunnelling), and is dominated by the capacitance operator. The other makes the transition between states of different 
wells, 8 Q — > o ± 7T, and only depends on the hopping irrelevant operator. This difference would help explain why the 
capacitance operator does not show up in the static conductance: its contribution to the oscillator strength is through 
transitions conserving O (charge). 
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In summary, we have studied the low energy properties of a LL in the presence of 2k f scattering in the repulsive 
regime g < 1. The numerical analysis of the spectrum for very large sizes has been made possible thanks to the 
use of the SCHA as an intermediate step that effectively integrates high energy states. Accuracy has been tested in 
the nominally worst case of small u. Presence of irrelevant operators and their manifestation in transport are two 
issues treated independently. We confirm the existence of two irrelevant operators, with g = 1/2 as the boundary 
between their regions of dominance. One represents particle hopping, and has dimension 1/g. The other describes 
charge oscillations around the pinned value, and has dimension 2. Temperature dependence of the static conductance 
has been show to be controlled by the hopping operator alone, giving G ~ T 2/>9 ~ 2 for all values of g. Frecuency 
dependence is affected by both operators, giving G ~ uj n , with rj = (2/g — 2) for g > 1/2, and ij = 2 for g < 1/2. 
Conceptual simplicity, tested numerical accuracy, and versatility in handling different aspects make this approach a 
promising tool for other quantum impurity problems. 

The author is grateful to Prof. F. Guinea for very stimulating discussions. This work has been supported by the 
DGICyT of Spain (grant PB92-0169). 
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FIG. 1. Two lowest bands £1,2(9) of the scaled hamiltonian h(L) = HL/tv with g — 0.4 for four increasing values of the 
scaled perturbation [Eq.(2)] u = + ,0.72, 1.63 and 34.9 (left to right). The continuous line of left panel is the exact result of a 
perfect LL. 



FIG. 2. Left panel: ln(8e\) [Eq.(3)] versus ln(L/ir) for g — 0.7,0.55,1/2,0.4, and 1/3 (top to bottom). Continuous lines 
are asymptotic behavior [Eq.(3)]. Right panel: as in left panel for ln(8e2) [Eq.(4)] with g = 0.7, 0.55, 1/3, 0.4, and 1/2 (top to 
bottom). 



FIG. 3. Left panel: ln(G) [Eq.(5)] versus ln(u) [Eq. (2)] for g = 0.4 with results for 200 (u,L) points collapsed onto a single 
scaling curve. Continuous line is the asymptotic behavior ~ u~ 2 ^ 9 . Arrow marks the exact limit for u = 0. Right panel: Zn(/i2) 
[Eq.(7)] versus ln(L/n) for g = 0.7, 0.55, 1/2, 1/3, and 0.4 (top to bottom). Continuous lines are asymptotic behavior [Eqs.(7)]. 
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